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Abstract 

A novel pedagogical technique is presented that can be used in the undergraduate (UG) class 
to formulate a relativistically extended Kinetic Theory of Gases and Maxwell-Boltzmann thermal 
speed distribution, while keeping the basic thermal symmetry arguments intact. The adopted 
framework can be used by students to understand the physics in a thermally governed system at 
high temperature and speeds , without having to indulge in high level tensor based mathematics. 1 



to/) 



Or 



Our approach will first recapitulate what is taught and known in the UG class and then present 
a methodology that will help students to understand and derive the physics of relativistic thermal 
systems. The methodology uses simple tools well known in the UG class and involves a compo- 
nent of computational techniques that can be used to involve students in this exercise. We also 
present towards the end the interesting implications of the relativistically extended distribution 
and compare it with Maxwell-Boltzmann results at various temperatures. 
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I. INTRODUCTION 



A. Thermal Gas in the Undergraduate Class 

The standard Kinetic Theory of Gases (KTG) as taught in a basic undergraduate (UG) 
course is a collision-based description of a classical ideal gas of particles, 2 with the walls of 
the container using a Newtonian framework and assumptions. 3 It results in expressions for 
the macroscopic properties of the gas (pressure, average kinetic energy K, root mean square 
speed f rms etc.) as a function of the absolute temperature T of the gas. The derivation of 
these results can be found in commonly used texts, 4 and yields the following major results for 
a classical gas of N particles each having a mass m with ks being the Boltzmann constant: 

v = — • (1) 
K = ^L. (2 ) 

A distribution function describes the probability of a particle's speed near a particular 
value as a function of the gas' absolute temperature, the mass of the particle and the value 
of speed under consideration. Random thermal motion dictates that both position and 
velocity space are uniformly distributed and the distribution function is stationary with 
time. A probability distribution function(PDF) f(.) of the speed of the particles in one 
direction must also give the directional probability in other independent directions as well. 
In the cartesian coordinate system, the probability of the speed of a particle to lie between 
Vk — > Vk + dvk for k = {x, y,z} is given by f(v x ).f(v y ).f(v z )dv x dv y dv z . Since any direction 
is as good as the others, the resultant distribution function must only depend on the 
total speed v of the particle, 

j(v 2 ) = M + vl + vl) = f(v x ).f(v y ).f(v z ) . (3) 

The form of the directional function f(.) is clearly that of an exponential as seen by tak- 
ing partial derivatives of Eq. (3) with respect to Vk, representing a Standard Probability 
Distribution function of the form, 

f(v k ) = Ae-** , (4) 

where A and B are constants depending on T and m. The Maxwell-Boltzmann Thermal 
speed distribution 5 (MB distribution) gives the probability P(v) dv of the particle speed to 



lie within an elementary volume dv x dv y dv z , or to lie between the speeds v and v + dv in 
the velocity space, centered on (v x ,v y ,v z ) , where dv x dv y dv z = 4n v 2 dv: 

P(v)dv = f(v x ).f(v y ).f(v z )dv x dv y dv z = 4irA 3 v 2 e~ Bv2 dv . (5) 

The constants A and B in Eq. (5) are determined using the two integration conditions: 

1. The Classical All Particle Condition 

Classically and non-relativistically, the speed v of all particles follows < v < oo, 
hence the PDF P(v) of Eq. (5) can be normalized as: 

pOO POD 

/ P(v)dv= / 4nA 3 v 2 e- Bv2 dv = 1 . (6) 
Jo Jo 

2. The Classical Kinetic Theory of Gases v rms Result 

Relating Eq. (1) with the probability distribution function of Eq. (5), we get, 

r°° r°° o ikT 

v *= / v *p( v )& v= / 4TrA 3 v A e- Bv dv = — . (7) 

Jo Jo rn 

Thus, 

- 4 = W = W ■ < 8) 

B. Need and Motivation to extend the Maxwell-Boltzmann Distribution 

The following points on the formulation of the MB distribution yield discrepancies in 
regard to the Special Theory of Relativity: 

I. The classical KTG is based on Newtonian mechanics and does not consider the effect 
of Special Relativity. The Normal probability distribution that the MB distribution 
uses in Eq. (5) is a gaussian and the variable, the speed of the particle v, extends till 
infinity. Classical MB distribution assumes that v e [ 0, oo ) as in the normalization 
scheme of Eq. (6), but Einstein's Special Theory of Relativity 6 limits v to c: the speed 
of light in free space. MB distribution at any T will always predict a non-zero fraction 
of particles to have speeds greater than c, which is physically forbidden and leads to 
inaccuracies in the distribution, particularly at high temperatures. 



2. As the particle gas approaches high temperatures, more particles will have v — > c and 
hence relativistic effects cannot be ignored. Thus, we cannot use K = (mv 2 /2) as is 
implied by Eq (1). Rather, using Special Theory, K — (7 — l)mc 2 where 

— 1 /2 

7 = (f — v 2 /c 2 ) is the Lorentz factor. The basic result of the KTG in Eq. (1) can 
be extended under the constraints of the Special Theory of Relativity as done in the 
next section. 
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FIG. 1. MB PDF (Solid Curve) and CDF (Broken Curve) for a classical electron gas at 
T = 10 9 K 

Astrophysical systems, e.g. the Intra cluster medium of galaxy clusters, can reach tem- 
peratures as high as 10 8 K, 7 which involve relativistic speeds and describing processes like 
Thermal Bremsstrahlung Emission 8 using MB distributions will lead to errors as depicted 
with the help of Fig. 1. It depicts the MB distribution for a classical electron gas at T = 10 9 K 
with v m ost probable = 0.58 c and shows that the Cumulative Distribution Function (CDF) ^ 1 
as v — > c and predicts about 3.2 % particles to have speeds beyond c. The error in the MB 
distribution is beyond negligence and is not an accurate description of the system. MB dis- 
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tributions are used extensively in the UG curriculum to describe different thermal systems, 
whose speeds can be relativistic as shown above and to allow for deeper understanding of the 
system at student level, we are motivated to extend the classical KTG and MB distribution 
into the relativistic regime by incorporating the postulates of Special Theory of Relativity. 

II. PEDAGOGICAL EXTENSION OF KTG INTO THE RELATIVISTIC REGIME 

We have achieved a novel approach to extend the KTG into the relativistic regime using 
Four Vectors 9 , a standard tool of the undergraduate physics class. This approach allows 
us to understand the physics behind the physical origin of the gas' macroscopic properties 
on similar lines of the standard KTG derivation using a collision based approach, under 
the constraints of Special Theory of Relativity. Previous works by Synge and Juttner, 110 
have developed a relativistic thermal speed distribution with the usage of tensor calculus, a 
mathematical tool not available with all UG students. Our methodology using Four Vectors 
uses minimal and elementary mathematics but offers insight into the physics of the problem. 

From an inertial frame of reference S, we observe a box of volume V = I 3 containing N 
classical particles with all standard KTG assumptions valid except that relativistic effects 
can be ignored. Let the i th particle in the box be described by the position four vector X = 
(x°, x, y, z). Let At be the proper time between events as measured in the particle's rest 
frame S[ and At be the improper time between events as measured in the frame S which 
are related by the Lorentz factor as: 

* = =->■ <») 

The corresponding acceleration four vector can be written in a form, 11 

A = (77C, 7 2 « + 77^) , (10) 

where v is the spatial velocity of the particle, a — dv/dt and 7 = 7 3 (a • v)/c 2 . 

We first model the particle-wall elastic collision in one dimension in the frame S and 
then generalize the results to three dimensions. In a one dimensional analysis along an axis 
perpendicular to a particular wall of the container, say x-axis: v = v x x and the acceleration 
a — dv/dt = Av/At. A similar approach as evolved in the Standard Newtonian KTG can 
be extended here such that on each collision with the particular wall of the container, we 
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have the following: 

Av = 2v x x, (11) 

At = - , (12) 

v x 

v 2 v 2 

a= u —x = -fx, (13) 

V V 

V 3 

a-v = -f. (14) 

Thus the spatial component of the acceleration four vector takes the form of, 
A _ {v 2 )ill) , K)(7, 4 ) _ {v 2 ){ll) u , fo 2 )^ 

Aspatial.i.lD - j 1 j-^ - j I 1 " 1 ~2 J' 

where ^ x — (1 — v 2 /c 2 ) is the Lorentz factor assigned in one dimension and 
since [1 + {v 2 {^l)/c 2 )] = (7^), hence we get: 

A _ / lfi x 
Aspatial,i,lD — ^ • UOJ 

The spatial part of the force four vector can be written as -F spa tiai,i,iD = m Aspatiai,i,iD and 
macroscopic gas properties can be derived as follows. The pressure due to the i th particle is 

nm- — y — • [i-i) 

The total pressure of the gas can be found by summing P^m in Eq. (17) over all the N 
particles in the box 

p to tai = e = e ^ mm , (is) 
1=1 i=i 

where (3 = v/c and dimension-specific beta (3k = Vk/c for k = {x, y, z}. 

The particles of the gas are in constant random motion under thermal equilibrium which 
provides that all directions away from the walls of the container are identical in all physical 
respects and no particular direction is preferred over the other. The symmetry in all the 
three cartesian directions due to thermal equilibrium of the gas leads to the conclusion that 
the average contribution of each cartesian direction component is the same, 

PI = P'y = ft = f , (19) 

where (3 X , f} y , f3 z are related to the beta factor j3 of the particle as, (3 2 = (3 2 + (3 2 + /3 2 . 

Mathematical analysis using infinite geometric series helps alternatively express Eq. (18) by 

introducing average of powers of j3 using Eq. (19) as: 

3 Mm r 2 

P^= 6 -^(P 2 W), (20) 

6 



where 



(21) 



As with the classical KTG, we assume that the gas particles behave ideally and follow the 
Ideal Gas Equation (PV = Nk B T) to generate the final Relativistic KTG postulate 



(22) 

Equation (22) is the relativistically extended version of the classical KTG postulate as in 
Eq. (1) by incorporating constraints of Special Relativity and can supersede the classical 
KTG postulate. 



III. FORMULATING A RELATIVISTIC DISTRIBUTION FUNCTION 

The underlying symmetry arguments in the MB distribution of uniform occupation of 
position and velocity space and equivalence among different directions as stated in §1 (A) 
are valid and are the defining criteria for random motion in thermal equilibrium. We can 
formulate the relativistically extended (RE) distribution function by incorporating changes 
in Eq. (5) compliant with Special Relativity. Let J{v)dv be the RE distribution function that 
has all the features of a thermal distribution. Postulates of Special Relativity will dictate 
that J(0) = J(v > c) = because no particle can travel with the vacuum speed of light or 
beyond. Eq. (3) will hold for the functional form of J{v)dv and the directional distribution 
function /re(-) can be assumed on similar lines of §1 (A) to be f-Rsivk) = Pe^ Rv ^ with P 
and R as constants. 

The resultant RE distribution function in velocity space centered on (v x ,v y ,v z ) or equiv- 
alently in (3 space centered on (j3 x , (3 y , (3 Z ) is, respectively, 

J(v)dv = AnP*v 2 e { - Rv ^ 2) dv , (23) 

J(f3)d(3 = 4nX 3 (3 2 e { - Y ^ 2) d(3 , (24) 

and the constants in these equations being related as X = Pc and Y = Rc 2 

It is important to mention at this stage that J(/3)dj3 is not strictly a Normal Probability 
Distribution 12 because the argument of the exponential in Eq. (24) does not vary as /3 2 but 
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as /3 2 7 2 i.e. (/3 2 /l — /3 2 ). Hence, Eq. (24) is a Probability Representative Function(PRF) 
and not strictly a PDF since it predicts a number proportional to the probability. The 
constants X and Y in Eq. (24) are dependant on T and m, and can be determined using the 
following two conditions: 

1. The Relativistic All Particle Condition 

As constrained by the Special Theory of Relativity, the speed v of any particle follows 
< v < c, hence the PRF J(f3)d(3 of Eq. (24) can be normalized to include all the N 
gas particles as: 

f J(f3)df3= [ 4nX 3 (3 2 e- Y ^ 2 d(3 = l. (25) 
Jo Jo 

2. The Relativistically Extended KTG Result 

Relating Eq. (22) and (21) with the PRF of Eq. (24), one gets: 

The integral conditions of Eq. (25) and (26) are not solvable in closed analytic form. The 
computational facility Wolfram Alpha 13 available was used to numerically solve these for a 
classical ideal gas of electrons by feeding a 40-value long array of constant Y to generate 
values of corresponding X and T. The Curve Fitting Tool of MATLAB 14 can be used to plot 
the temperature dependance of X and Y (Fig. 2) and apply a best fit for the data wherever 
possible. A linear polynomial best fit the correlation of ln(X) - ln(T) and the T dependance 
of Xis 

In X= 10.4- (0.5063) In T. (27) 



IV. FEATURES OF THE RELATIVISTICALLY EXTENDED DISTRIBUTION 
A. The Most Probable Speed Analysis 

The most probable (mp) speed of the particles in the MB distribution is v mp = 
(2/cT/m) 1 / 2 . A similar expression using the RE distribution can be found when the PRF of 
Eq. (24) is maximum, which gives 

(i4i P ) 2 = F' (28) 
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FIG. 2. Temperature dependance of constants X (Solid Curve) and Y (Broken Curve) of Eq. (24) 

which can be solved using a quadratic analysis of roots. In the non-relativistic limit v « c 
so that 7 — y 1, 2 « (1/V) = (l/i?c 2 ) from Eq. (24) and (23), which is equivalent to 
v^ p w (1/-R). The constant i? in the MB PDF of Eq. (5) is ^ which gives us, 

2A;T 

Vp,v«c = ( ) > (29) 
m 

and hence in the non-relativistic regime, the results for most probable speed using the RE 
distribution match the MB results. 

B. The nature of temperature dependance of constant Y 

The argument of the exponential in the PRF of Eq. (23) must be unit-less, hence R 
will have the units of (sec 2 /m 2 ) or equivalently (kg/joule) i.e inverse of specific work. R 
in Eq. (23) is inversely proportional to the work or energy required to heat the gas to the 
temperature T. The temperature dependance of Y as in Fig. 2 is analogous to the 
x-axis reflected plot of total energy E = (mc 2 / yT— /3 2 ) of a particle with it's {3 parameter. 
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As the temperature and the required energy of the gas to accomplish this increases there is 
decrease in Y. After a certain T range, any further increase in T requires an uncontrollable 
increase in energy input and a corresponding uncontrollable decrease in Y. It is clear that 
a high T involves relativistic speeds and nature demands lot of energy to sustain such high 
speeds and temperatures. 



V. A COMPARISON STUDY 



We now compare results of the MB and RE distributions at various temperatures and 
show that results of MB distribution for T > 5 x 10 8 K are inaccurate beyond negligence and 
the RE distribution results are physically convincing and can be a good approximation to 
the real system. As stated in §111, Eq. (24) is a PRF and predicts a number proportional to 
the probability. The MB and RE distributions can be compared at a given T by a parameter 
(r) at various speeds from Eq. (5) and Eq. (23). 

Ratio(r) = ^4 • (30) 
J(y) 

We will use r max to represent the Ratio (r) at the highest plotted speed and r m i n will represent 
the Ratio (r) at the lowest plotted speed and Ar = r max — r m j n , an average estimate of the 
scatter between the two functions at a given T. 
Table I chalks out a comparison (Fig. 3 and 4) between the MB and RE distribution 

TABLE I. Comparing MB and RE distribution functions at two different temperatures 



T = 2.96 x 10 5 K 


T = 9.18 x 10 s K 


r max = 3.52 x KT 11 


rmax = 4.31 X 10 9 


r min = 3.53 x KT 11 


r min = 4.16 x 10- 11 


Ar = 1.86 x KT 12 


Ar = 4.31 x 10 9 


^(mp-MB) = 0.0099984c 


f(mp-MB) = 0.5565418c 


*Vp-re) = 0.009999c 


^(mp-RE) = 0.618034c 


100 x f™P(v) dv = 


100 x f™P(v)dv = 3.2% 


J c °°J(v) dv = 


J c °°J(v)dv = 



functions at two different temperatures by comparing values of various parameters, 
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FIG. 3. MB PDF (Solid Curve) and Relativistically Extended (RE) PRF (Broken Curve) for 
T = 2.96 xlO 5 K 



representative of each distribution. As evident from the comparison table, for 
T^5x 10 8 K, MB and RE distributions match very well with the shape of the distribution 
function preserved and and the scatter parameter r between the two is of the order of 1CT 12 
or less. RE distribution converges to Maxwellian results at low speeds and low temperatures. 

Figures 3 and 4 plot MB and RE distribution functions at the two comparison tempera- 
tures. They clearly show that for higher T, the shape of the two functions varies significantly, 
with MB predicting about 3.2% particles beyond c, while the RE function is a good descrip- 
tion of the system under the constraints of special relativity. Figure 5 shows the variation 
of the Ratio rwith j3 for T = 2.96 x 10 5 K. The scatter in the two distributions is extremely 
large as v — > c but it is important to notice that at all temperature regimes, RE results 
converge to MB at low speeds. Overall, the RE and MB distributions diverge at high T and 
high v and the RE distribution favors higher speeds compared to MB, as shown by the v mp . 
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FIG. 4. MB PDF (Solid Curve) and Relativistically Extended (RE) PRF (Broken Curve) for 
T = 9.18 xlO 8 K 

VI. CONCLUSION 

In our analysis of including Special Relativity in the KTG and MB distribution, we have 
presented a pedagogically useful method to introduce relativistic constraints that preserve 
the basic collision and thermal framework that is taught as the standard approach to this 
subject. The method of Four Vectors and particle-wall collisions along with the usage of 
computational techniques to solve integral conditions can be very effectively used in the 
UG classroom paradigm to help students appreciate the interlink between different areas of 
Physics and themselves try and formulate it. The ideas presented here give a reasonable 
description of the physical system but it is important to note that this approach is based on 
incorporation of Special Relativity in a previously accepted methodology and the formulation 
of the RE function involves assuming a logically sound function, extrapolating from the MB 
framework. Thus, we suggest that due care must be taken while adopting our approach 
which we feel is a good pedagogical tool in the UG classroom. 
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